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A powerful mathematical method for the investigation of the properties of dynamical 
systems is represented by the Kosambi-Cartan-Chern (KCC) theory. In this approach 
the time evolution of a dynamical system is described in geometric terms, treating 
the solution curves of a dynamical system by geometrical methods inspired by the 
geodesics theory of Finsler spaces. In order to geometrize the dynamical evolution 
one introduces a non-linear and a Berwald type connection, respectively, and thus 
the properties of any dynamical system are described in terms of five geometrical 
invariants, with the second one giving the Jacobi stability of the system. The Jacobi 
(in)stability is a natural generalization of the (in)stability of the geodesic flow on a 
differentiable manifold endowed with a metric (Riemannian or Finslerian) to the non¬ 
metric setting. Usually, the KCC theory is formulated by reducing the dynamical 
evolution equations to a set of second order differential equations. In this paper we 
present an alternative view on the KCC theory, in which the theory is applied to a 
first order dynamical system. After introducing the general framework of the KCC 
theory, we investigate in detail the properties of the two dimensional autonomous 
dynamical systems. The relationship between the linear stability and the Jacobi 
stability is also established. As a physical application of the formalism we consider 
the geometrization of Hamiltonian systems with one degree of freedom, and their 
stability properties. 
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I. INTRODUCTION 


One of the fundamental theoretical tools extensively used in every branch of natural 
sciences for modeling the evolution of natural phenomena are the continuously time evolving 
dynamical systems. In scientific applications they are extremely useful, and their usefulness 
is determined by their predictive power. In turn, this predictive power is mostly determined 
by the stability of their solutions. In a realistic setup some uncertainty in the measured 
initial conditions in a physical system always does exist. Therefore a mathematical model 
that is physically meaningful must offer information and control on the time evolution of 
the deviations of the trajectories of the dynamical system with respect to a given reference 
trajectory. It is important to note that a local understanding of the stability is equally 
important as the global understanding of the late-time deviations. From a mathematical 
point of view the global stability of the solutions of the dynamical systems is described by the 
well known theory of the Lyapunov stability. In this approach to stability the basic quantities 
are the Lyapunov exponents, measuring exponential deviations from the given reference 
trajectoryii^. However, one should mention that it is usually very difficult to determine the 
Lyapunov exponents analytically. Therefore, various numerical methods for their calculation 
have been developed, and are applied in the study of the dynamical systems^-d^. 

However, even though the methods of the Lyapunov stability analysis are well established 
and well understood, it is important to adopt in the study the stability of the dynamical 
system different points of view. Then, once such a study is done, one can compare the 
obtained alternative results with the corresponding Lyapunov linear stability analysis. An 
important alternative approach to the study of the dynamical systems is represented by 
what we may call the geometro-dynamical approach. An example of such an approach 
is Kosambi-Cartan-Chern (KCC) theory, which was initiated in the pioneering works of 
Kosambi^, Cartanl^ and Cherni^, respectively. The KCC theory is inspired, and based, 
on the geometry of the Finsler spaces. Its basic idea is the fundamental assumption that 
there is a one to one correspondence between a second order dynamical system, and the 
geodesic equations in an associated Finsler space (for a recent review of the KCC theory 
see^). From a geometric point of view the KCC theory is a differential geometric theory of 
the variational equations describing the deviations of the whole trajectory of a dynamical 
system with respect to the nearby ones^^. In this geometrical description to each dynamical 
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system one associates a non-linear connection, and a Berwald type connection. With the 
help of these two connections hve geometrical invariants are obtained. The most important 
of them is the second invariant, also called the curvature deviation tensor. Its importance 
relies on the fact that it gives the Jacobi stability of the systen>i^^— . The KCC theory 
has been extensively applied for the study of different physical, biochemical or engineering 
systemsiSi^— . 

An alternative geometrization method for dynamical systems was introduced irt^ and^^, 
and further developed in~2I. Applications to the Henon-Heiles system and Bianchi type 
IX cosmological models were also investigated. In particular, iii^ a theoretical approach 
describing geometrically the behavior of dynamical systems, and of their chaotic properties 
was considered. In this case for the underground manifold a Finsler space was adopted. The 
properties of the Finsler space allow the description of a wide class of dynamical systems, 
including those with potentials depending on time and velocities. These are systems for 
which the Riemannian geometry approach is generally unsuitable. 

The Riemannian geometric approach to dynamical systems is based on the well-known 
results that the flow associated with a time dependent Hamiltonian 

+ ( 1 ) 


can be described as a geodesic flow in a curved, but conformally flat, manifold^. With 
the introduction of a metric of the form = W (a;“) Sabdx°'dx^, in which the conformal 
factor is given by W {x°') = E — V {x°') , where E is the conserved energy associated with 
the time-independent Hamiltonian if, it follows that in the metric Qab = W (a;“) dab-, the 
geodesic equation for motion is completely equivalent to the Hamilton equations^ 


dx°' OH dpa 


dH 


( 2 ) 


dt dpa' dt dXa 
This result implies that the confluence, or divergence of nearby trajectories x“(s) and [x + 
.^]“(s) of the Hamiltonian dynamical system is determined by the Jacobi equation, i.e., the 
equation of geodesic deviation, which for the present case takes the following form. 


= RUuVe = -KX- 




( 3 ) 


In Eq. ([3]) is the Riemann tensor associated with the metric Qab, and D/Ds = u°‘Va 
denotes a directional derivative along the velocity held = dx'^/ds. Therefore linear 
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stability for the trajectory x“(s) is thus related to the Riemann curvature or, more 
exactly, to the curvature K^. If, for example, is everywhere negative, then it follows 
that K'^ always has one or more negative eigenvalues, and therefore the trajectory must be 
linearly unstable^^. 

There are a large number of mathematical results on the geometrization of the dynamical 
systems. For example, from a geometric perspective point of view, iii^S the global and local 
stability of solutions of continuously evolving dynamical systems was reconsidered, and the 
local stability was defined based on the choice of a linear connection. Note that an important 
point in favor of the use of a linear connection is that it is naturally dehned for any dynamical 
system (S', X), and not only for those related to second-order evolution equations. 

An important testing ground of the KCC theory is represented by the study of two- 
dimensional autonomous systems, and of their stability properties. Such a study was per¬ 
formed for two dimensional systems of the form 

under the assumption that the point (0,0) is a hxed point, i.e. /(0,0) = 5'(0,0) = 0. By 
relabeling v as x, and g{u,v) as y, and by assuming that 5'u|(o,o) 7^ 0, one can eliminate the 
variable u. Moreover, since {u,v) = (0,0) is a fixed point, from the Theorem of Implicit 
Functions it follows that the equation g{u, x) —y = 0 has a solution u = u{x, y) in the vicinity 
of (x, y) = (0, 0). Since x = g = gu f+gv 2/, we obtain an autonomous one-dimensional second 
order equation, equivalent to the system (|1]), namely 

+ g^{x,y) = 0, (5) 

where 

g\x, y) = -gu{u{x, y),x) /(m(x, 2 /), x) - gMx, y),x) y. (6) 

The Jacobi stability properties of Eq. ([5]) can be studied by using the KCC theoryJ^J^, and 
the comparative study of the Jacobi and Lyapunov stability can be performed in detail. 

In the present paper we will introduce an alternative view of the KCC theory by adopting 
the perspective of the hrst order dynamical systems. After a brief presentation of the gen¬ 
eral formalism, we will concentrate our attention to the simple, but important case of two 
dimensional autonomous dynamical systems, whose properties are studied in detail. Instead 
of reducing the two-dimensional autonomous system to a single second order differential 


5 


equation of the form ([5]), by taking the time derivative of each equation we transform the 
system to an equivalent system of two second-order differential equations. We also clarify 


the relation of the Jacobi stability approach with classical Lyapunov stability theory. As a 
physical application of our formalism we apply it to the study of two dimensional Hamil¬ 


tonian systems, describing physical processes with one degree of freedom. The Hamilton 
system is transformed into an equivalent system of two second order differential equations, 
which can be studied geometrically similarly as geodesic equations into an associated Finsler 
space. We obtain the geometrical quantities describing the geometerized Hamilton system, 
and the conditions for its Jacobi stability are obtained. 

The present paper is organized as follows. The KCC theory is presented in Section [ITl 
The applications of the KCC stability theory to the case of two-dimensional systems is 
considered in Section m The comparison between the linear Lyapounov stability and the 
Jacobi stability is performed in Section |IVl The application of the KCC theory to two 
dimensional Hamiltonian systems is discussed in Section |Vl We discuss and conclude our 
results in Section ED 

II. BRIEF REVIEW OF THE KCC THEORY AND OF THE JACOBI 
STABILITY 

In the present Section we briefly summarize the basic concepts and results of the KCC 
theory, and we introduce the relevant notations (for a detailed presentation see^^ ancJl^). 

A. Geometrical interpretation of dynamical systems 

In the following we assume that AJ is a real, smooth n-dimensional manifold, and we 
denote by TAJ its tangent bundle. On an open connected subset H of the Euclidian (2?7,-|-1) 
dimensional space i?” x x R^, we introduce a 2n -|- 1 dimensional coordinates system 
{x\ ?/*, t), i = 1, 2,..., n, where (A) = ..., x”), (j/*) = y ‘^,..., y'^) and t is the usual 

time coordinate. The coordinates r/® are dehned as 



(7) 


A basic assumption in our approach is the assumption that the time coordinate t is an 
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absolute invariant. Therefore, the only admissible coordinate transformations are 


i = x”) ,i G {1,2, 


(8) 


Following2^, we assume that a deterministic dynamical systems can be dehned as a set of 
formal rules that describe the evolution of points in a set S with respect to an external time 
parameter t ^ T, which can be discrete, or continuous. More exactly, a dynamical system 
is a map^ 

(j) : T X S ^ S, {t, x) ^ x), (9) 


which satishes the condition 0(t, •) o 0(s, •) = (j){t + s, •), Vt, s ^ T. For realistic dynamical 
systems that can model natural phenomena additional structures need to be added to the 
above dehnition. 

In many situations of physical interest the equations of motion of a dynamical system 
follow from a Lagrangian L via the Euler-Lagrange equations. 


d dL dL 
dt dy^ dx^ 


1,2, ...,n. 


( 10 ) 


where Fi, i = l,2,...,n, is the external force. Note that the triplet {M,L,Fi) is called a 
Finslerian mechanical systen>^>^. If the Lagrangian L is regular, it follows that the Euler- 
Lagrange equations dehned in Eq. flTOj) are equivalent to a system of second-order ordinary 
(usually nonlinear) differential equations 


dF 


+ 2G^ {x^ 


y\t) = 0,z e (1,2, ...,n} , 


( 11 ) 


where each function G'^{x^,y^,t) is G°° in a neighborhood of some initial conditions 
iix)oAy)oFo) in fl. 

The fundamental idea of the KCC theory is that if an arbitrary system of second-order 
differential equations of the form ffTTj) is given, with no a priori Lagrangean function as¬ 
sumed, still one can study the behavior of its trajectories by analogy with the trajectories 
of the Euler-Lagrange system. 


B. The non-linear connection and the KCC invariants associated to a 
dynamical system 

As a hrst step in the analysis of the geometry associated to the dynamical system dehned 
by Eqs. (na, we introduce a nonlinear connection N on M, with coefficients Aj, dehned 
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(12) 




dG^ 


m = 

^ dyi 

Geometrically the nonlinear connection iVj can be interpreted in terms of a dynamical 
covariant derivative for two vector fields n, w defined over a manifold M, we define the 
covariant derivative as^ 


V% = 


+ N]ix, y)w^ 
dx^ ' 


A 

dx^ 


(13) 


For Nl{x, y) = Tlj^{x)y\ from Eq. (11311 we recover the dehnition of the covariant derivative 
for the special case of a standard linear connection, as dehned in Riemmannian geometry. 

For the non-singular coordinate transformations introduced through Eqs. (jS]), the KCC- 
covariant differential of a vector held on the open subset G C x x is dehned 




dt dt ^ 


For = y^ we obtain 


Dy^ 


= Ny - 2G* = -eh 


(14) 


(16) 


dt 

The contravariant vector held e* dehned on G is called the hrst KCC invariant. 

Now we vary the trajectories of the system (ITTl) into nearby ones according to the 
rule 

R (t) = x\t) RyCit), (16) 

where |? 7 | is a small parameter, and ^*(t) are the components of a contravariant vector held 
dehned along the trajectory x’‘(t). By substituting Eqs. (ITHll into Eqs. (ITTll . and by taking 
the limit r/ —)■ 0, we obtain the deviation, or Jacobi, equations in the formr^^— 


dR 




( 17 ) 


Eq. (flTll can be rewritten in a covariant form with the use of the KCC-covariant derivative 


as 


where we have denoted 


dR 


= py, 


p; = -2 


dG^ , ■ ,dm . , dm 


(18) 


(19) 
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(20) 


and we have introduced the Berwald connection G*;, dehned as^— 



dN] 
dy’' ' 


The tensor Pj is called the second KCC-invariant, or the deviation curvature tensor, 
while Eq. flTSD is called the Jacobi equation. When the system of equations (ITT]) describes 
the geodesic equations, in either Riemann or Finsler geometry Eq. ffTSjl is the Jacobi held 
equation. 

The trace P of the curvature deviation tensor can be obtained from the relation 


P = P! = 


- 2G'G‘, + + N;nI + 


( 21 ) 


dx’- ' ‘ 

One can also introduce the third, fourth and hfth invariants of the system ffTTl) . which 
are dehned asl^ 


P* = - 

jk 3 


1 zap* dPt 


TP ^ 

’ ^jkl r. I 1 ^jkl 


dyk Qyi^-3ki Qyi ■ ( 22 ) 

Geometrically, the third invariant P^f, can be interpreted as a torsion tensor. The fourth 
and hfth invariants Pjj^i and are called the Riemann-Christohel curvature tensor, and 
the Douglas tensor, respectivelyl^J^. Note that in a Berwald space these tensors always 
exist. In the KCC theory they describe the geometrical properties and interpretation of a 
system of second-order diherential equations. 


C. The definition of the Jacobi stability 

The behavior of the trajectories of the dynamical system given by Eqs. (ITTil in a vicinity of 
a point X® (to) is extremely important in many physical, chemical or biological applications. 
For simplicity in the following we take to = 0. In the following we consider the trajectories 
X® = x®(t) as curves in the Euclidean space (P®®, (.,.)), where (.,.) is the canonical inner 
product of P®®. We assume that the deviation vector Z obeys the initial conditions Z (0) = O 
and Z (0) = JE 7 ^ O, where O G P®® is the null vector—^—. 

Thus, for the focusing tendency of the trajectories around to = 0 we introduce the fol¬ 
lowing description: if ||Z(^)|| < ~ O’*", then the trajectories are bunching together. 

But if ||Z(^)|| > t ~ O’*", the trajectories have a dispersing behavior—^—. The focus¬ 

ing/dispersing tendency of the trajectories of a dynamical system can be also described in 
terms of the deviation curvature tensor in the following way: The trajectories of the system 
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of equations (fTT]) are bunching together for t O’*" if and only if the real part of the eigenval¬ 
ues of the deviation tensor Pj (0) are strictly negative. On the other hand, the trajectories 
are dispersing if and only if the real part of the eigenvalues of Pj ( 0 ) are strictly positivei^'— . 

Based on the above considerations we introduce the concept of the Jacobi stability for a 
dynamical system as follows^^— : 

Definition: If the system of differential equations Eqs. (1TT|1 satishes the initial conditions 
||P (to) — P (0)11 = 0 , ||P (to) — P (0)11 7^ 0 , with respect to the norm ||.|| induced by a 
positive dehnite inner product, then the trajectories of Eqs. (fTTj) are Jacobi stable if and only 
if the real parts of the eigenvalues of the deviation tensor Pj are strictly negative everywhere. 
Otherwise, the trajectories are Jacobi unstable. 


III. THE CASE OF THE TWO VARIABLE DEPENDENT DYNAMICAL 
SYSTEM 


In the present Section we consider the case of the two-dimensional dynamical systems 
already studied inJ^J^. However, the present formulation is in a parametric form. 

Let’s consider the following two dimensional dynamical system, 

(23) 

^ = g [x\x^) . (24) 

We point out that we regard a solution of Eqs. fl25]l and fl24D as a flow 93^ : P C —)■ 
or, more generally, 93 * : P C M —)■ M, where M is a smooth surface in The canonical 
lift of (ft to the tangent space TM can be geometrically dehned as 


(ft '■ TM ^ TM, (25) 

ipt{u) = {ipt{u),(pt{u)). (26) 

In terms of dynamical systems, we simply take the derivative of Eqs. (|23|) and fl24)) and 
obtain 


(fx^ 

dt'^ 

d‘^x‘^ 


/l (x\ 1/^ +/2 (x\ 1/^ 

(27) 

9i + g2 {x^,x^) y^, 

(28) 
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where we have denoted 


dx^ 2 


(29) 


t — t ®9 . ®9 .,1 _ 

* O 1 ’ ^* 0 2’ * O 1 ’ * 02’^ ’ y ^4- 

ox^ ox^ ox^ ox^ dt dt 

In other words, on TM we obtain 

d^ 
dt 

dt 

where {x^x‘^,y^,y‘^) are local coordinates on TM. Hence we can see that the above system 
is actually a linear dynamical system on the hber T(a,i 3 ; 2 )M. Moreover, the system of Eqs. 
(I27D and fl28|l can be written as 

d^x^ 


fi {x\x‘^) y^ + /s {x\x‘^) y^, 

(30) 

91 {x^,x^) y^ + 92 {x\x‘^) 2/^ 

(31) 


df^ 

d^x"^ 

dt"^ 


+ (-/i2/^ - f2y‘^) = 0, 

+ - W) = 0. 


By comparison with Eqs. fllip we have 


(32) 

(33) 



ilhv' + hy^ 




where 


J = J{f.g) = 


is the Jacobian of the dynamical system Eqs. 
nonlinear connection 



(34) 


(35) 


. Using Eq. flT^ we obtain for the 


1 US)-(11) -5 («»)■+"■” « 

Therefore all the components of the Berwald connection cancel, 

dm 

g;, ^ (37) 

Then, for the components of the deviation curvature tensor (Ej), given in Eq. ([19]), we 


obtain 


{dG^\ \dNl odN} ,.1 ^ ^ r 1 


dx^ 2^ dx^ ' ' 


= ^W + ^/i22/' + iV/iVj, 
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= -2 


f dG‘^\ idN? Ar2AH ^ / 1 2\ 1 1 2 ^ 5 '! Ar2AH 

(&rj+» ^+!' ^ (sis +» )- 2 !/ N,, 


and so on. Therefore we get 


{11) = 


pi pi 
-'1 -'2 


p2 p2 
1 -'2 

/iii/^ + fi2y^ giiy^ + gi2y^ 
fi2y^ + f22y^ gi2y^ + 9229 "^ 


1 / fiiy^ + /i 22 /^ fuy^ + / 22 i/^ 


+ fi'l 2 l/^ 5'12|/^ + 9229^ 

t 


+ -^JIU.9)^J\U\9) 


+ \jlU.9)><J\U.9) 


fii /i2\ (y^ 
J2I f22 W 


9ii 912 \ y 


1 


+ ^Jiif,g)xJlif,g). 


^ 5 ' 2 i g 22 J \y 

Therefore we have obtained the following 

Proposition 3.1 The curvature deviation tensor associated to a second order dynamical 
system is given by 


P = 


{Hf V H,-y)‘ + \r-U, 9 ), 


(38) 


where Hf = 


fll fl2 
/ 2 I / 22 , 


is the Hessian of f, and similarly for g. 


IV. LYAPUNOV AND JACOBI STABILITY OF TWO DIMENSIONAL 
DYNAMICAL SYSTEMS 


Similarly within, and withont losing generality, we assnme p = (0, 0) is a hxed point of 
Eqs. fl2^ and fl2T|) . Then the Lyapnnov stability is governed by the characteristic eqnation 

/i^ — tr A ■ /i + det A = 0, (39) 


where tr A and det A are the trace and determinant of the matrix 


^ Jify g)\{o,o) 



(40) 


We also denote by A = (tr A)^ — 4 det A the discriminant of Eq. fl3^ . By descending 
again on (or M), and evalnate at the hxed point (0,0), we obtain |(o,o) = (0,0), 

and therefore 


^ 1 ( 0 , 0 ) = 



(41) 
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To be more precise, we have 


fl /2\ f fl /2\ if /l + /25'1 / 1/2 + 1292 
9i 92 j \gi 92 J ^ yfi9i + 9i 92 f29i + 92 

Hence we have 

H P = - (/i + 2/291 + 92 ) = {fl + 2 / 1(72 + 92 ~ 2 / 1(72 + 2 / 291 ) = - [(tr Af — 2 det A) , 

(43) 

det P = ^ (/i ^2 - / 29 if = Q det . (44) 

The eigenvalues of the matrix (P*) are given by the characteristic equation 




— tr P ■ A + det P = 0, 


(45) 


and its discriminant is 

A = (tr P)^—4detP = — [(tr A^ — 2 det Hi ^ — ^ (detH)^ = — (tr H)^ [(tr H)^ — 4 det Hi . 

16 4 16 

(46) 

Thus we obtain the following 

Computational Lemma 4.1. The trace, determinant and discriminant 0 /the charac¬ 
teristic equation 0 / the deviation curvature matrix P are: 

tr P = ^ [(tr H)^ — 2 det H] , 

det P = — (det Af , 

16 

A = ^(tr A)'A, 

where A = (tr H)^ — 4 det H is the discriminant 0 / Eq. 

We recall now some general results of linear algebra. 

Lemma 4.2. Let H be a (2,2) matrix, and denote by Ai, A 2 its eigenvalues. Then 

(i) the eigenvalues of fc • H are k ■ \i, k ■ A 2 , for fc 7 ^ 0 scalar. 

(ii) the eigenvalues of H^ := H • A.. • H are (Ai)*^, (A 2 )^. 

From here it follows: 

Lemma 4.3. 1/ Xi, A 2 are eigenvalues 0 / A, then pi = (|Ai)^, 7 x 2 = (|-^ 2 )^ are eigen¬ 
values 0 / P. 
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Remark: The formulas of Lemma 4.3 imply 


S' /ii + p2 — ^ (Ai + A 2 ) — - [(Ai + A2)^ — 2 ■ A 1 A 2 ] , 


p Pi ■ /i 2 — ( -Ai 


H =S(W 


The above formulas are consistent with the computational Lemma 4.1. 


A. Comparison of Jacobi and Lyapunov stability 

Definition Let p be a fixed point of the two-dimensional system x = /(x), and denote 
by Ai, A 2 the two eigenvalues of A := (Df)ip. The following classihcation of the hxed point 
p is standard. 

(I) A > 0, Ai, A 2 are real and distinct. 

1 . 1 . Ai ■ A 2 > 0 (the eigenvalues have the same sign): p is called a node or type 
I singularity; that is, every orbit tends to the origin in a dehnite direction as 

t —)■ CXD. 

1.1.1. Ai, A 2 > 0 : p is an unstable node. 

1 . 1 . 2 . Ai, A 2 < 0 : p is a stable node. 

1.2. A < 0, Ai • A 2 < 0 (the eigenvalues have different signs): p is an unstable 
fixed point, or a saddle point singularity. 

(II) Ai, A 2 are complex, i.e. Ai ,2 = o ± ifi, fi 0- 

11 . 1 . a 7 ^ 0 : p is a spiral, or a focus, that is, the solutions approach the origin as 
t —)■ cxD, but not from a dehnite direction. 

11 . 1 . 1 . a < 0 : p is a stable focus. 

11.1.2. q; > 0: p is an unstable focus. 

11 . 2 . 0 = 0 : p is a center, that means it is not stable in the usual sense, and we have 
to look at higher order derivatives. 

(Ill) A = 0, Ai, A 2 are equal, i.e. Ai = A 2 = A. 
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111.1. If there are two linearly independent eigenvectors, we have a star singularity, or 
a stable singular node (these are simple straight lines through the origin). 

111.2. If there is only one linearly independent eigenvector, we have an improper node, 
or unstable degenerate node. 

By combining the above results with the Computational Lemma 4.1, it follows that 

I. A > 0, /ii,/r 2 £ 

1.1. S > 0, p is Jacobi unstable 

(tr AY — 4 det A > 0, (tr — 2 det A > 0. (47) 

- for det A > 0, we must have (tr A)^ > 2 det A > 4 det A, and hence both cases 
(tr A > 0, det A > 0) and (tr A < 0, det A < 0) imply Jacobi stability. 

-for det A < 0, we obtain again Jacobi stability because now relations fl47|) are identically 
satished. 

1.2. S' > 0: this case is not possible algebraically. 

II. A < 0: /ii ^2 = a ± i/5 e C (Remark that A < 0 ^ A < 0, i.e. Ai, A 2 are complex, and 
hence det A > 0 . 

11.1. S > 0: p is Jacobi unstable 

(tr AY — 4 det A < 0 , (tr A)^ — 2 det A > 0 , (48) 

that is, 2 det A < (trA)^ < 4 det A. This case is possible for (tr A > 0,detA > 0), and 

(tr A < 0 , det A > 0 ). 

For det A > 0 this case is not possible algebraically. 

11.2. S < 0: p is Jacobi stable 

(tr A)^ — 4 det A < 0 , (tr A)^ — 2 det A < 0 , (49) 

that is (tr A)^ < 2 det A < 4 det A. This is possible for both (tr A > 0 ,detA > 0), and 

(tr A < 0, det A > 0). 

III. A = 0: pi = P 2 £ equivalent to A = 0, or tr A = 0. 

- if A = 0, then Ai = A 2 € M and p is singular node or degenerate node 
-if tr A = 0, A 7 ^ 0, it follows 

(i) A > 0, det A < 0, i.e. saddle point, or 

(ii) A < 0 , det A < 0 , i.e., p is a center. 
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B. Important remarks 


1. Assume p is Jacobi stable, that is, by definition, we must have one of the following 
situations: 

(i) A > 0, S' > 0, or 

(ii) A < 0, S > 0. 

As we have seen already (i) is algebraically impossible, hence, if p is Jacobi stable it must 
follow A < 0 -v^ A < 0. Hence we have proved: 

If p is Jacobi stable, then A < 0. 

2. Conversely, assume A < 0. This is equivalent to A < 0 due to the Computational 
Lemma 3.3. Since A < 0 implies det A > 0, and Ai ^2 = a ± i/9, we have 

/^i = (^Ai) =^{a + il3f = ^[{a^ - 13'^)+2iaj3], (50) 

h2=QA2^ = ^[{a^ - -2iaP]. (51) 

From the above equations we obtain 

s = 1 (a= - . (52) 

Therefore it follows 

(i) — /9^ > 0 ^ p is Jacobi unstable, 

(ii) — /9^ < 0 ^ p is Jacobi stable. 

Condition (ii) above does not appear inJ^ due to the reduction of the two dimensional 
dynamical system to a one-dimensional SODE. We obtain 

Theorem 4.1. 1. If p is a Jacobi stable fixed point, then A < 0. 

2 . If A < 0 for the fixed point p, then 

(i) if — /9^ > 0, then p is Jacobi unstable; 

(ii) if — /9^ < 0, then p is Jacobi stable, where Ai _2 = a ± i/9 are the eigenvalues of A. 

V. APPLICATIONS TO TWO DIMENSIONAL HAMILTONIAN SYSTEMS 

As a physical application of the formalism developed in this paper, in the present Section 
we consider the geometrical description, and the stability properties of a physical system 
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described by a Hamiltonian function H = H{x,p) = H where H{x,p) G C"", n > 2. 

From a physical point of view x represents the particle’s coordinate, while p is its momentum. 
The motion of the system is described by a Hamiltonian system of equations in the plane. 

Definition. A system of differential equations on is called a conservative Hamiltonian 
system with one degree of freedom if it can be expressed in the form 


dx^ dH{x^,x‘^) 


dt 


dx‘^ 


= H 2 {x^,x'^) , 


dx"^ dH{x^,x‘^) 


dt 


dx^ 




(53) 


By taking the total derivative of the Hamiltonian function, and with the use of Eqs. 
we obtain 

dH dH {x^,x'^) dx^ dH{x^,x'^)dx‘^_ 
dt dx^ dt dx"^ dt 

Eqs. fIMl) shows that H {x^,x^) is constant along the solution curves of Eqs. fl53|) . Hence 


(54) 


the Hamiltonian function is a hrst integral and a constant of motion. Moreover, all the 
trajectories of the dynamical system lie on the contours dehned by H {x^,x‘^) = C, where C 
is a constant. From a physical point of view H {x^,x‘^) represents the total energy, which is 
a conserved quantity. 

By taking the derivatives of Eqs. fl53|) with respect to the time parameter t we obtain 
hrst 

d'^x^ 
dt"^ 

d'^x'^ 
dt"^ 


: H 21 {x^,x‘^) + H 22 {x^,x‘^) 

(55) 

-Hn y^ - Hi 2 (a;\ i/^. 

(56) 


which can be written in the equivalent form 

d^x 


dt‘^ 


+ 2G' = 0,! = 1,2, 


(57) 


where 


/gM 

l=-l( 

Igu 

1 2( 



(58) 


where 


Jh = Jh{-Hi,H2) = 


H. 


21 


H. 


22 


-Hu -H 


12 , 


(59) 
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is the Jacobian of the Hamiltonian system given by Eqs. fl5^ . Eqs. flFTI) give the geometrical 
interpretation of a bidimensional (one degree of freedom) Hamiltonian system, showing that 
they can be studied by similar methods as the geodesics in a Finsler space. 

Using Eq. flT^ we obtain for the nonlinear connection associated to a Hamiltonian system 
the expressions 




Nl 

iVf A^l 



= (60) 


Therefore, for a Hamiltonian system all the components of the Berwald connection vanish. 


dm 


(61) 


The components of the deviation curvature tensor of a Hamiltonian system can be ob¬ 
tained as 



H 


(\ ^ ( miiy^ + H 2 W + H222y^ 

p 2 j ^ y-Hniy^ - P1121/' -P1122/' - Hu 2 y^ 

^iJH):{-m,H2)xiJH)]i-m,H2) 



t 

+ 


+ 


1 

4 


(Jh)] (-Pi, P 2 ) X (Jh)] (-Pi, P 2 ) • 


Therefore we have the following 

Proposition 5.1 The curvature deviation tensor associated to a Hamiltonian dynamical 
system is given by 


Ph = 1 (h.„, -y + \jl (--ffi, , 


(62) 


P2II P2I2 

where = is the Hessian of Hi, and similarly for H 2 - 

-Hill -Pii2^ 

Hence in the following we can introduce the concept of the Jacobi stability of a Hamilto¬ 
nian system by means of the following 
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Definition. If the Hamiltonian system of Eqs. 


satisfies the initial conditions 


||x* (to) — 5^* (^o)ll = 0, I |x* (to) — (^o) 11 7^ 0, with respect to the norm ||.|| induced by 

a positive definite inner product, then the trajectories of Hamiltonian dynamical system are 
Jacobi stable if and only if the real parts of the eigenvalues of the curvature deviation tensor 
Ph are strictly negative everywhere. Otherwise, the trajectories are Jacobi unstable. 

To illustrate the implications of the geometric approach introduced here we consider the 
simple case of the one dimensional conservative motion of a point particle with mass m > 0 
under the influence of an external potential V{x), described by the Hamiltonian 

.2 1 


p- 


= — + V(x) = — (xY + H (x^) . 


(63) 


The equation of motion of the particle are 


dX^ 1 n dx'^ / IX 

-IP = ’ -PP = ’ 

dt m dt 


(64) 


where in the following we denote by a prime the derivative with respect to the coordinate 
x^. By taking the derivative with respect to time of Eqs. fIM)) we obtain first 


d'^x^ 1 dx'^ 


d^x^ 


dP 


m dt ' dP 




(65) 


Eqs. 


can be written as 


d‘^x^ 

dP 


+ 2G^ 2/^) = 0, i = 1, 2, 


( 66 ) 


where 


G‘ {x\x^ y\ v^) = -Y’ !,=) = \v" (x‘) y'. 


The components of the non-linear connection iVj = dG^/dy^ can be obtained as 
.u dG^ r. dG^ 1 ..2 dG^ Ixv/. n ..2 dG^ ^ 

' ~ Thp ~ ' ~ ‘ “ 2 '"^ h ) ■ 2 - - ■ 


(67) 


( 68 ) 


2m dy^ 2 

With the use of Eq. (IT^ we obtain for the components of the deviation curvature tensor 


P* the expressions 


, Pi = 0, Pi = --V" (x^) 2/\ P| = -^V" (xM 

1 4m ^ ^ ^ ^ 2 ^ ^ ^ ’ 2 4m ^ ^ 


(69) 


The eigenvalues Ai 2 of the curvature deviation tensor are given by 


Ai ,2 - 2 


Pi + P| ± ^APlPl + {Pi - Pif 


(70) 
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and can be obtained explicitly as 


A.. - - p.). 


(71) 


Therefore we have obtained the following 

Theorem 5.1 The trajectories of a one-dimensional Hamiltonian dynamical system with 
point particle like Hamiltonian given by Eq. fl6^ are Jacobi stable if and only if the potential 
V{x) of the external forces satisfies for all x the condition V'\x) > 0. 

From a physical point of view the condition V"{x) > 0 implies that the potential V{x) 
is in a minimnm. More exactly, if V (xq) = 0 is an eqnilibrinm state of a physical system 
at point Xo, the condition for the Jacobi stability V" (xq) > 0 reqnires that the potential V 
has a minimnm at x = xq. 

For a point particle like two dimensional Hamiltonian system with one degree of freedom 
the eqnations of the geodesic deviation, satisfied by the deviation vector i = 1,2, given 
by Eq. (im. can be obtained as 


1 

dt'^ m dt 


d^e 

dt'^ 


+ V" (x^) 


dt 


+ V" (i‘) y'e = 0. 


(72) 

(73) 


Let’s assnme now that (x^ = xo,x^ = 0,y^ = 0,y‘^ = 0) are the critical points of the 
Hamiltonian system of Eqs. fIMl) . Then the system of the geodesic deviation eqnations 
take the form 


^ 

dt'^ m dt ’ 


fe 

dt'^ 


+ V" (xo) 


dt 


0 , 


(74) 

(75) 


and mnst be integrated with the initial conditions ^^(0) = 0, ^^(0) = 0, = .^lo, and 

^^(0) = ^ 20 ) respectively. The geodesic deviation eqnations for the point mass Hamiltonian 
have the general solntion 


eit) 


1 

V" (xo) 


(^o) sin 


\dJMtj 

Vm j 


- 6o cos 


y/m j 


+ '^20 


(76) 


e{t) 




cos 


\/m j 



+ 


\/^ 6 o 


sin 


y/m J ' 
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(77) 
























VI. DISCUSSIONS AND FINAL REMARKS 


In the present paper we have considered an alternative approach to the standard KCC 
theory for hrst order autonomous dynamical systems, based on a different transformation 
of the system to second order differential equations. The approach was presented in detail 
for two dimensional dynamical systems, for which the two basic stability analysis methods - 
the (Lyapunov) linear stability analysis and the Jacobi stability analysis - were discussed in 
detail. From the point of view of the KCC theory the present approach allows an extension 
of the geometric framework for hrst order systems, such increasing the parameter space, and 
the predictive power, of the method. We have also found that there is a good correlation 
between the linear stability of the critical points, and the Jacobi stability of the same points, 
describing the robustness of the corresponding trajectory to a small perturbation^^. On the 
other hand, the Jacobi stability is a very convenient way of describing the resistance of limit 
cycles to small perturbation of trajectories. 

As an application of our approach we have considered the study of a bi-dimensional 
Hamiltonian system, describing the one dimensional (one degree of freedom) motion of a 
physical system. The KCC theory can provide an alternative, and very powerful, method 
for the geometrization of classical mechanical systems, whose properties can be described 
by a Hamiltonian function. The transformation of the corresponding Hamilton equations 
to second order differential equations allows naturally their study similarly as geodesics in 
an associated Finsler space, and gives the possibility of a full geometric description of the 
properties of the dynamical system in a non-metric setting. This represents one of the basic 
differences, and advantages, of the KCC approach as compared to the alternative Jacobi 
and Eisenhart methods for geometrization^^, which essentially require as a starting point 
a metric. It is important to emphasize that the advantages of the geometric approach to 
the description of the dynamical systems are not only conceptual, but also the method has 
a predictive value. By starting from the deviation curvature tensor we can obtain some 
effective stability conditions for physical systems. Moreover, in the present approach, the 
geodesic deviation equation can be formulated and solved rather easily (either analytically or 
numerically), and thus the behavior of the full perturbations of the trajectories near critical 
points can be studied in detail. 

To summarize our results, in the present paper we have introduced and studied in de- 
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tail some geometrical theoretical tools necessary for an in depth analysis and description 
of the stability properties of dynamical systems that may play a fundamental role in our 
understanding of the evolution of natural phenomena. 
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